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ABSTRACT 

A grid of evolutionary sequences of stars in the mass range 1.2-7M0, with solar-like 
initial composition is presented. We focus on this mass range in order to estimate 
the masses and calculate the CNO surface abundances of a sample of observed red 
giants. The stellar models are calculated from the zero-age main sequence till the early 
asymptotic giant branch (AGB) phase. Stars of M ^ 2.2M0 are evolved through the 
core helium flash. In this work, an approach is adopted that improves the mass de¬ 
termination of an observed sample of 21 RGB and early AGB stars. This approach 
is based on comparing the observationally derived effective temperatures and abso¬ 
lute magnitudes with the calculated values based on our evolutionary tracks in the 
Hertzsprung-Russell diagram. A more reliable determination of the stellar masses is 
achieved by using evolutionary tracks extended to the range of observation. In ad¬ 
dition, the predicted GNO surface abundances are compared to the observationally 
inferred values in order to show how far standard evolutionary calculation can be 
used to interpret available observations and to illustrate the role of convective mix¬ 
ing. We find that extra mixing beyond the convective boundary determined by the 
Schwarzschild criterion is needed to explain the observational oxygen isotopic ratios 
in low mass stars. The effect of recent determinations of proton capture reactions and 
their uncertainties on the and ratios is also shown. It is found that 

the ^^N(p, 7)^®0 reaction is important for predicting the ratio in red giants. 

Key words: convection, nuclear reactions, nucleosynthesis, abundances - stars: evo¬ 
lution - stars: low-mass 


1 INTRODUCTION 

After the main sequence evolutionary phase, stars evolve 
to the red giant branch (RGB). This evolution is initiated 
by the ignition of shell H-burning surrounding the He core, 
whose energy flux causes the envelope to expand and the star 
evolves to the RGB. The expansion increases the opacity and 
leads to the development of a deep convective envelope. This 
is the first dredge up event (FDUP), as convection mixes up 
the products of H-burning to the surface altering the surface 
composition of the star. 

In the mass range (4-7)M0, stars exhibit blue 
loops at the beginning of core He burning (see 
Halabi. El Eid fc Champagne! (l2012l l. and references 
therein). The main phase of core He-burning is completed 
before the track evolves back to the RGB. This leads 
again to the deepening of envelope convection. For solar 
metallicity stars of M >4 Mq, a second dredge up (SDUP) 
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can reach deeper regions, which leads to further changes in 
the surface abundances. 

This wor k uses observ ations obtained for a sample of 
red giants by ITsu U (Hooi), hereafter TsujiOS, in order to 
achieve two goals: (a) to estimate the masses of these giants 
by matching their observationally derived effective tempera¬ 
tures and bolometric magnitudes to the values obtained from 
the evolutionary tracks in the HR diagram. This is possible 
since the stars are not pulsating Mira variables (see Section 
3.1 for details). We are able to improve the mass determina¬ 
tion done by TsujiOS by using more extended evolutionary 
tracks to avoid the extrapolation that he partially relied on 
to determine the mass of some giants, (b) to compare the cal¬ 
culated CNO abundances of these models to those inferred 
from observations in the light of recent determinations of 
key reaction rates. 


A large body of observational data is available for the 
surfa ce CNO abundances i n RGB st ars dLambert fc R.eisl 


198ll: Harris & Lambert 1984al bl: iHarris et al.l 

1988; 

Lambert et al.l Il98fil: ICilrov & BrownI Il99ll: iTsuiil 

1991 ; 
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Charbonnell 1 1994 iTsuiil l2008l : iTautvaisiene et aP I 2 OI 0 I : 
Piau et alj boill b These data provide a powerful tool to 


get insight into the internal structure of evolved stars. 
In particular, comparing the predicted oxygen surface 
abundances to the observed data is useful to improve the 
treatment of convective mixing in the stellar interiors (see 
Section 3.2.2). Our results suggest the need for extra mixing 
below the edge of the convective envelope as determined 
by the Schwarzschild criterion in order to achieve a better 
agreement with observations. This extra mixing has been 
suggested in several investi gations in connection with the 
evolu t ion of field g i ants ifCharbonnel fc Do Nasciment 


Tautvaisiene et al.l 


199^^ Gratton_at_^J_ 2000ll. open clusters | Lu^ 


n^td 

199i 


2003 

: Pilachowski et al. 

2003 

2007 

: iDenissenkov et al. 

2 OI 5 I) 


200oj, 2005 ), globular clustem dShetroM 


200 s ; IPecio-Blanco fc de Lavernv 
I 2 OI 5 II as well as to explain i s otopic 


IBusso et al.l l2014l l. The abundance profile of of par¬ 
ticular interest. This is because this isotope is produced 
by the ON-cycle which requires higher temperatures than 
the CN-cycle. Therefore, the profile exhibits a steep 
gradient within the central region of the star, at the position 
of maximum convective penetration (shown later in Fig. 
Hi. This renders the surface abundance sensitive to 
the depth of convective mixing, stellar mass and to the 
nucle ar reaction rates involved in the CNO cycle (IEI EidI 
1 1994) . We show in Section 3.2.2 how the ratio 

is useful to constrain this extra mixing. The low surface 
carbon isotopic ratios in low mass stars however, can 
not be explained by our extra-mixing treatment. Other 
non-standard mixing me chanisms may need to be inv oked 
as shown in Section 3.2. iRamstedt fc OlofssonI ll20l4 ) also 
discuss carbon isotopic ratios in AGB stars in connection 
with their recent estimations of circumstellar ^^CO/^^CO 
abundance ratios based on radiative transfer analysis of 
radio line emission observations. 

Concerning the nuclear reaction network, there has been 
extensive experimental work recently on updated determi¬ 
nations of major reaction rates, inclu ding those of the CNO 
cycle. One set of rates that we use llSergi et al.| 2014 ) was 
obtai ned with the Trojan Horse method ( La Cognata et ^ 
l20inl) . which is an indirect technique that is able to provide 
more reliable reaction rate cross-sections at low tempera¬ 
tures where measurements in the astro physical energy range 
are available. Other se ts, obtained bv ISallaska et aH (|2013l ) 
and llliad is et al. 1 I 2 OIC 1 I) are evalua ted based on Monte Carlo 
techniques i Loneland et al ] l2010l) . This method provides a 
median rate which -under certain conditions- resembles the 
commonly referred to “recommended” rate, as well as a low 
rate and a high rate which, unlike the “upper” and “lower” 
limits of classical rates, have a well-defined statistical mean¬ 
ing. The effect of these recent determinations of the proton 
capture reactions including the ^^N(p, 7)^®0 rate on the iso¬ 
topic ratios ^®0/^^0 and is investigated. 

This paper is organized as follows. Section 2 describes 
the calculations and provides a summary of the observa¬ 
tional data used in this work. The results are presented in 
Section 3. Mass determinations and comparison with pre¬ 
vious works are provided in Section 3.1. In Section 3.2, we 
describe the surface abundance profiles and the effects of ex¬ 
tra mixing. The effect of nuclear reaction rates is discussed 
in Section 3.3. Concluding remarks are given in Section 4. 


2 MODEL CALCULATIONS 


The evolutionary sequences presented in this work are ob¬ 
tained us ing the stellar evo lution code HYADES as de¬ 
scribed in |ELEid_et_alJ (l2()09|), with the re c ent m odifications 
outlined in Halabi. El Eid fc Champagnel ll2012l) . This code 
is a one-dimensional implicit Lagrangian code based on a 
hydrodynamical method which solves the stellar structure 
equations on an adaptive grid. Mass-loss is included using 
semi-empirical rates adjusted to the global parameters of the 
star. The rates are used according to a Mira puls ation pe¬ 
riod j P) llVassiliadis fc Wo^ll99^ . For P<100 d. iReime^ 
|l97a) mass-loss rate is used, with 7 = 1 . For 100 ^P<500 d 
we use a more effective rate according to iBowe 3 (Il988h . For 
P ^ 500 d, the super wind mass-loss rate during the ACB is 
used as suggested bv I Vassiliadis fc WoodI (Il992f ). 


2.1 Analyzing the overshooting region 


In the context of the Mixing Length Theory (MLT), 
the extension of a convective zone is determined by the 
Schwarzschild criterion, that is when where 

and Vgjj are the radiative and adiabatic temperature 
gradients, respectively. 

Within this framework, a long-standing issue is to fix 
the boundary of the convective zone. In the local description 
of the MLT, mixing beyond the Schwarzschild boundary is 
introduced in a parameterized way. Accordi ng the multi¬ 
dimen sional hydrodynamic simulations by iFrevtag et al.l 
lll996l) . a local description of an extra mixing (or overshoot¬ 
ing) may be introduced in terms of an exponentially decay¬ 
ing diffusion coefficient: 

-2Z 

D(z) = (1) 


where z= I’^boundary “ overshoot distance. Do 

is the diffusion coefficient at the boundary of the convec- 
tive envelope obtain ed from the Mixing Length Theory (see 
iLanger et al.l (Il985l )). Hp is the pressure scale height and / 
is a free parameter which is a measure of the efficiency of this 
extra partial mixing. It is clear from Eq. [T] that for smaller 
values of /, D has a steeper profile, or equivalently less ex¬ 
tra mixing. As / increases, this extra mixing extends further 
beyond the formal convect i ve bou ndary. The numerical sim¬ 
ulations by IFrevtag et al.l ill 9961 ) find / =0.25 ± 0.05 and 
1.0 ± 0.1 for A-stars and DA white dwarfs, respectively. We 
will use the observationally inferred oxygen isotopic ratio in 
red giants to constrain the value of / (see Section 3.2). 

The mixing of chemical elements is achieved by solving 
the diffusion equation given by: 


dX, 

dt 


d 

dMr 


[(47rr^p)^D 


dXi 

dMr 


( 2 ) 


where r is the radius, p is the mass density and D is the 
diffusion coefficient given by Eq. [T] when used in the over¬ 
shoot region, otherwise it is equal to Do in a convective zone 
according to the Schwarzschild criterion. 

To illustrate the effect of the treatment of mixing de¬ 
scribed above. Fig. [1] shows the behavior of Vj-gj , Vg^^ 
as well as the profiles of the diffusion coefficient and hydro¬ 
gen (Xjj) as a function of interior mass. This is done for a 
1.2Mq model during the first dredge up phase (FDUP) after 
the star has evolved to the red giant branch. 
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Figure 1. Physical quantities related to the condition of convection at the bottom of the convective envelope in a 1.2M0 star during 
FDUP, for (a) the standard stellar model and (b) the model with overshooting using f=0.125. (c) and (d) are the same as (a) and (b), 
respectively, but 3x10^ years later. and are the radiative and adiabatic temperature gradients, respectively. Also shown are 

the hydrogen profile (left scale) and diffusion coefficient profile (right scale). 


The left-hand panel of Fig. [T] shows the extension of 
the convective envelope in the standard models, i.e. without 
overshooting where the abrupt drop of the diffusion coeffi¬ 
cient D is visible along with the vertical drop of Xjj. The 
right panel shows the behavior with overshooting. It is seen 
that in the latter case, mixing is extended into the radia¬ 
tive region where The lower panels show the 

profiles at the next time-step, 3x10® years later. In both 
cases, the convective envelope deepens in mass as FDUP 
proceeds, but it is deeper in the case with overshooting, 
where the change in composition alters the opacity so that 
^rad ^ ^ad becomes satished deeper in mass. Therefore, 
overshooting does not only induce extra mixing but also 
drives convective instability. Details of the present calcula¬ 
tions with this overshooting or extra diffusive mixing are 
provided in Section 3.2.2. 


2.2 Observational data 

A sample of red giant stars has been observed by TsujiOS as 
given in Table [T] The effective temperatures were ob tained 
using the infrared flux method jBlackwell et al.lll98di . while 
the absolute bolometric magnitudes were determined from 
the bolometric luminosities obtained by integrating the spec¬ 
tral energy distributions and the Hipparcos parallaxes. The 
uncertainty on the effective temperature is estimated to be 


lOOK, and the error on the bolometric magnitude is mainly 
due to the error on the parallaxes. In the next section, the 
effective temperatures and bolometric magnitudes are used 
to determine the stellar masses of the observed giants us¬ 
ing our evolutionary tracks, which cover the whole range of 
observations. 


3 EVOLUTIONARY RESULTS 


Evolutionary sequences of stars in the mass range (1.2- 
7 . 3 )M 0 are calculated from the zero-age main sequence till 
the early AGB phase. Stars of M ^ 3Mq exhibit blue loops 
starting at the onset of central helium burning, which be¬ 
come more extended for stars of mass > dM©. A detailed 
di scussion on this evolutionary phase has been presented 
in [Halabi. El Eid fc Champagnel (l2012l l. Stars of masses ^ 
2 . 2 M 0 and sol ar-like initial composition evolve through the 
core He-flash ||Kippenhahn fc Weiger^ llQQOl : iMocak et al.l 
I 2 OIOI : iBildsten et alJ 2012li . We hnd that for M < 2 M 0 , the 
helium flash starts off-center owing to the cooling via the 
plasma and photo neutrino energy losses. In the mass range 
2 ^M/M 0 2.2 the helium flash starts at the center since 

these stars evolve at relatively lower central densities so the 
cooling by neutrino energy losses is less efficient. 

The core helium flash requires very short time steps to 
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Table 1. Spectral type, effective temperatures and bolometric 
magnitudes (TsujiOS). 


Object(BS/HD) 

Spectral type 

Te// (K) 

Mbol (mag) 

5 Vir (4910) 

M3HI 

3643 

-2.4 ± 0.3 

a. Tau (1457) 

K5-I-HI 

3874 

-1.7± 0.2 

RRUMi (5589) 

M4.5HI 

3397 

-3.4 ± 0.3 

RZ Ari (687) 

M6IH 

3341 

-3.5 ± 0.6 

5 Oph (6056) 

M0.5HI 

3790 

-2.2 ± 0.3 

n Vir (4517) 

Mini 

3812 

-2.2 ± 0.4 

Eri (1003) 

M3-l-HIa 

3712 

-2.9 ± 0.4 

10 Dra (5226) 

M3.5HI 

3730 

-2.9 ± 0.3 

/3 Peg (8775) 

M2.5H-HI 

3606 

-3.3 ± 0.2 

305 Her (6146) 

M6-HI 

3298 

-4.2 ± 0.4 

(T Lib (5603) 

M2.5HI 

3596 

-3.4 ± 0.5 

R Lyr (7157) 

M5HI 

3313 

-4.3 ± 0.3 

g Gem (2286) 

M3HI 

3643 

-3.3 ± 0.3 

OP Her (6702) 

M5H 

3325 

-4.4 ± 0.8 

p Per (921) 

M4H 

3523 

-4.1 ± 0.4 

a Get (911) 

M1.5HIa 

3909 

-3.2 ± 0.3 

A Aqr (8698) 

M2.5HI 

3852 

-3.4 ± 0.7 

XY Lyr (7009) 

M5H 

3300 

-5.1 ± 1.1 

5^ Lyr (7139) 

M4H 

3420 

-5.5 ± 0.8 

a Her (6406) 

M5Ib-H 

3293 

-5.8 ± 1.6 

BS6861(6861) 

M4 

3600 

-5.2 ± 2.0 


accommodate the rapidly changing variables. In our calcula¬ 
tion, the time step is of the order of less than a year during 
the core helium flash. For stars of masses in excess of 2 . 2 M 0 , 
no significant degeneracy effects occur and core He-burning 
proceeds under hydrostatic conditions. 


3.1 Mass evaluation 


In the work by TsujiOS, the masses of the red giant stars 
listed in Tabl e [J are derived using the evolutionary tracks by 
IClared ll2004h. A main difference betwee n our code and that 
of [a ared ll2004h is that the latter used ICauehlan fc Fowler] 
lllOSSh rates for the basic nuclear reactions in the network, 
while our used reaction rates are updated accord ing to 
the JINA REACLIB database llCvburt et al.|[2010l l. This 


is expected to introduce modifications on the evolutionary 
tracks. Another difference is that Claret assumes core over¬ 
shooting, with an overshooting distance of 0.2Hp for the 
whole mass range and ignores envelope overshooting. Core 
overshooting results in a bigger core mass and affects the 
evolutionary tracks and the stellar lifetimes. In the sam¬ 
ple considered here, most of the stars are low-mass stars 
which have very small or no convective core at all, thus ap¬ 
plying the same amount of overshooting at the convective 


core boundary would result in a large amount of mixing 
that yields results which are inconsistent with observations 
dWoo fc Demarau^ l200lh . In order to avoid this artifact, 
core overshooting is not incl uded in our c alculation. 

Moreover, the tracks bv I Clare t| ||2004|) are not extended 
enough, so that extrapolations were needed at temperatures 
below 3200K during the RGB phase for about half of the 
stars in the sample studied by TsujiOS. In this calculation, 
the stars are evolved until the early AGB phase without 
relying on extrapolation. Moreover, mass-loss is included, 
and we report the masses of the evolved stars rather than 
their initial masses, which improves the mass determination, 
particularly for the higher masses where mass-loss becomes 
more effective. 

Adopting the values of Mi^^j and Te// given in Ta¬ 
ble [T] the evolutionary tracks shown in Fig. [5] are used to 
evaluate the stellar masses. It is important to note here 
that a direct comparison of the theoretical temperature 
to that inferred observationally wouldn’t have been possi¬ 
ble if the stars are pulsating Mira variables, in which case 
a radius cannot be strictly defined and any comparison 
wouldn’t hold. Even the term effective temperature may be¬ 
come questionable for the very evolved AGB stars featuring 
strong pulsations an d mass loss rates llBaschek et al.lll99ll : 
iLebzelter et al.ll2010ll . However, in this sample, the stars are 
on the RGB or early AGB phase, and thus, haven’t yet expe¬ 
rienced any thermal pulsations. This allows a reliable com¬ 
parison between the predicted effective temperatures and 
the observationally inferred ones. The theoretical bolomet¬ 
ric magnitude is obtained using the well known relation: 
Mboi = 4.75-2.5log(L/L 0 ). 

As seen in Eig. |3] our evolutionary tracks describe well 
the advanced evolutionary stage of these stars. Having ob¬ 
tained the mass of every star using these tracks, it is possible 
to identify its evolutionary stage and compare its GNO sur¬ 
face abundances with the observational data. This will be 
described in Section 3.2. 

The stellar masses are given in Table [2] Eor complete¬ 
ness, m asses obtained by other works for some stars are in¬ 
cluded (lMaillardlll974l:IS mith fc Lamberjl 19851 : [Harris et al.l 
1 19881: iDecin et al.lll997ll . The error on the mass is deter¬ 
mined from the error bars on the observational M^^^j and 
Te//. Table [2] shows that our values are systematically lower 
than those by TsujiOS. We attribute this mainly to two rea¬ 
sons: 

(a) The stars are evolved to the stage where they are 
observed, that is, we do not use any extrapolated tracks as 
done in TsujiOS. 

(b) Mass-loss is taken into consideration, which becomes 
significant for the more massive stars. 

Moreover, it is clear from Table [2] that our errors on the 
masses are also lower. Calculating the tracks up to advanced 
stages helps to get better evaluation of the masses of red 
giants. 


3.2 CNO Surface Abundances in Red Giants 

3.2.1 Predictions of surface abundances with standard 
mixing 

This section presents the results for the surface CNO abun¬ 
dances of the studied sample of stars. These are then com- 
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Figure 2. Bolometric magnitude versus effective temperature (Te//) showing the present evolutionary tracks of stars of masses 

(1.2-7)M0. Also shown is the observed sample listed in Tabled 



log TeffiK) 


Figure 3. The evolutionary tracks together with the observed red giants listed in Tabled Note that the observed data points have 
errors as indicated in Tabled 
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Table 2. a: present work is the mass on the RGB, in¬ 

cluding mass-loss), b: TsujiOS, c: Smith and Lambert (1985), d: 
Harris et al. (1988), e: Decin et al. (1997), /: El Eid (1994). The 
masses are given in solar units. Note that the error on the mass 
by Tsuji08 increases for more massive stars, which may be related 
to the inaccuracy in the parallax measurements. 


Object(BS/HD) 

M “ 

initial 

M “ 

M 

M(others) 

5 Vir (4910) 

1.2 ± 0.2 

1.19 ± 0.2 

1.4 ±0.3 

2.0 = 

a Tau (1457) 

1.2 ± 0.2 

1.20 ±0.2 

1.5 ±0.3 

1.5 

RRUMi (5589) 

1.3 ± 0.1 

1.15 ±0.1 

1.6 ±0.3 


RZ Ari (687) 

1.3 ± 0.2 

1.14 ±0.2 

1.5 ±0.4 


5 Oph (6056) 

1.4 ± 0.2 

1.39 ±0.2 

1.6 ±0.3 


1 / Vir (4517) 

1.4 ± 0.4 

1.39 ±0.4 

1.7 ±0.4 

2.0 = 

Eri (1003) 

1.8 ± 0.3 

1.73 ±0.3 

2.0 ±0.4 


10 Dra (5226) 

1.8 ± 0.3 

1.74 ±0.3 

2.1 ± 1.8 


/3 Peg (8775) 

1.8 ± 0.3 

1.70 ±0.3 

2.2 ±0.3 

1.7 cA.eJ 

30g Her (6146) 

1.8 ± 0.3 

1.65 ±0.3 

2.0 ±0.6 

4.0= 

(T Lib (5603) 

2.0 ± 0.2 

1.90 ±0.2 

2.2 ±0.5 


R Lyr (7157) 

2.0 ± 0.2 

1.80 ±0.2 

2.1 ±0.5 


fj. Gem (2286) 

2.0 ± 0.5 

1.90 ±0.5 

2.3 ±0.5 

2.0 

OP Her (6702) 

2.1 ± 0.4 

1.90 ±0.4 

2.3 ± 1.0 


p Per (921) 

2.5 ± 0.4 

2.40 ±0.5 

3.2 ±0.5 


a Get (911) 

3.0 ± 0.5 

2.96 ±0.5 

3.6 ±0.4 


A Aqr (8698) 

3.0 ± 0.5 

2.96 ±0.5 

3.7 ± 1.2 


XY Lyr (7009) 

4.0 ± 1.0 

3.0 ± 1.0 

3.7 ± 1.5 


52 Lyr (7139) 

5.0 ± 1.0 

4.5 ± 1.0 

5.5 ± 2.0 


a Her (6406) 

5.5 ± 1.5 

4.0 ± 1.0 

5.0 ±2.0 

7.0 f 

BS 6861(6861) 

5.5 ± 1.5 

5.35 ± 1.5 

6.3 ±4.0 



pared to the CNO isotopic ratios inferred from observations. 
The abundance profiles of the CNO isotopes and those of H 
and He prior to the FDUP are shown in Fig. |4] In order to 
study the variation of these abundance profiles as a func¬ 
tion of the initial stellar mass, two sample stellar masses are 
shown: l.dM© (Fig.|4^) and SMq (Fig.|4jD). Several features 
can be identified in these figures: 

(a) The isotope is produced near the middle of the 
star by the CN cycle in all cases. This reflects the relatively 
low temperatures required for the production of this isotope 
by the CN-cycle. 

(b) The profile is remarkable, showing a steep gra¬ 
dient in the central region. This is because is produced 
by the ON cycle which requires higher temperatures to be¬ 
come effective. 

(c) The isotope is fragile and effectively destroyed 
by the ^®0(p, q)^®N reaction. 

These results are well known in the literature, but it 
is important to understand the surface abundances result¬ 
ing after FDUP and SDUP for different stellar masses. The 


efficiency of FDUP in modifying the surface abundances is 
related to the maximum penetration of the convective enve¬ 
lope on the RGB. In Fig. 01 the solid vertical line marks the 
convective boundary as determined by the Schwarzschild cri¬ 
terion and the dashed vertical line marks the position of the 
boundary when extra mixing beyond the formal convective 
boundary is considered (see Section 3.2.2 for details). 

After FDUP, the change in the surface composition of 
a certain isotope depends on the shape of its profile inside 
the star. For example, the peak of the profile is located 
near the middle part of the star, so that envelope convection 
is able to smear out the profile, causing an increase in the 

surface abundance, or a decrease in the ratio. 

For the isotope the situation is highly dependent on 
the stellar mass because the main production of is con¬ 
centrated in the inner part. Fig. |4^ which shows the case of 
a 1.4M0 as an example, illustrates that in low mass stars 
(M^ 2 M 0 ) envelope convection does not completely smear 
out the peak. This makes the profile sensitive to 
mixing in stars of 2 M 0 , so that any additional mixing 
below the envelope will increase its surface abundance. This 
effect is less pronounced in stars with M 3 M 0 (Fig. |4}3 
shows the case of a 5 M 0 as an example in this mass range), 
where the bump is fully engulfed by the formal convec¬ 
tive envelope, so extra mixing will not significantly alter its 
surface abundance in this mass range. 

Table |3] summarizes the values of and 

i'‘N/i®N after the FDUP and SDUP (if any, as we show 
later) in the case of standard convective mixing, together 
with those inferred from observations by TsujiOS and other 
independent field star observations. The stars in the sample 
of TsujiOS are advanced in evolution, but did not experience 
the third dredge up during the AGB phase. This is evident 
from the carbon surface abundances (Tsuji 2014, private 
communication). Thus, the comparison can be restricted to 
the effects of FDUP and SDUP only. 

It is not easy to make a direct comparison between 
predicted surface abundances and observations, since this 
comparison is model-dependent on theoretical and observa¬ 
tional grounds. Deriving the isotopic ratios observationally 
involves several sources of uncertainty like the dispersion in 
the ratios obtained from different lines, and inaccuracy in 
the atmospheric model parameters. Systematic errors may 
also be present, such as the uncertainty in the continuum po¬ 
sition a nd departures fro m local thermodynamic equilibrium 
(LTE) (lAbia et al.ll2^) . in addition to the difficulties that 
are inherent in the spectral analysis of very cool stars. The 
fact that atmospheric values are model-dependent causes 
discrepancies between observational results among different 
groups, and consequently, a ffect the subsequent discussion 
l|R,amstedt fc OlofssonI 120141 ). On the other hand, theoreti¬ 
cal models are also challenged by uncertainties in convective 
mixing, mass loss and nuclear reaction rates, where standard 
FDUP models often face difficulties in explaining carbon and 
oxygen surface abundances, particularly in low-mass stars. 
While being vigilant to these limitations, a careful compar¬ 
ison is useful for the sake of a better understanding. 

Fig. [5] shows the ratios as a function 

of stellar mass, along with the th eoretical pr ed iction s 
by iBoo throvd fc Sackm annl (|l999h . lAbia et ^ (|2ni2l l. 
iKarakas fc Lattanzid (l2014l ) and the available observations. 
The predicted shows a distinct behavior between 
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Mr/Ms Mr^s 


Figure 4. Abundance profiles of H, He and several CNO elements after core H-burning in the stellar masses (a) 1.4M0 and (b) 5 Mq. The 
solid lines mark the deepest penetration in mass of the convective envelope, while the dashed lines mark that when envelope overshooting 
is invoked. See text for details. 


high and low mass stars, which will be discussed separately. 
For the higher masses (M> 2 M 0 ), the standard calculation 
predicts ratios that ht the observed range within 

the error bars on the stellar masses. The relative spread in 
the observational values in this mass range cannot be ex¬ 
plained by overshooting since in these stars the bump in 
the *^^0 proHle is fully mixed by the formal convective enve¬ 
lope as shown in Fig. |4 Jd, so overshooting will not introduce 
a signihcant change to the standard predictions. Therefore, 
we attribute this spread to the uncertainties in describing 
the properties of the relatively cool red giants. This is a 
general feature of observationally inferred data due to the 
above mentioned uncertainties in measuring faint lines in 
their spectra which may result in considerably variant ratios 
for the same object among different observations. We expect 
that the error bars of these ratios are underestimated, par¬ 
ticularly in stars like a Get and a Her whose ratio varies 
considerably among different works as shown in Table O 

In the case of the lower mass stars, the theoretical 
isq/i^o ig higher than most data points obtained from 
observations. The results of standard calculation by other 
groups shown in the figure indicate a similar behavior. It is 
noted that the high observational values in the low 

mass stars like a Tau, v Vir, cr Lib, P Peg, 5 Vir, and RR 
UMi, which are denoted by up-pointing triangles in Fig. [5l 
are estimations due to the uncertainties in measuring faint 
lines like in relatively cool atmospheres, where at 

such low temperatures the absorption lines are strong caus¬ 
ing severe blending by several weak lines and introducing un¬ 
certainties (TsujiOS). Due to such complications, weak lines 
couldn’t be measured at all for in these stars, and 

isq/i^o ig not well determined. Thus, these observational 
data points are not reliable to compare to our calculation. 
Our discussion is based on the low observational ratios in low 
mass stars which are more dehnitive, with well-determined 
error bars. 

The overproduction of in low mass stars in our 

model shows that more needs to be mixed to the surface 
in order to lower the ratio. This may be achieved by 


present work: FDUP 
present work: SDUP 
Boothroyd&Sackmann(1999) 
Abiaetal.(2012) 
Karakas&Lattanzio(2014) 
TsujiOS 

Tsu]i08(lower limit estimates) 
Harris&Lambert(1984b) 
Harris&Lambertjl 988) 
Harris&Lambert(1985) 
Maillard(1974) for a Her 
Abia+(2012) for a Tau 



2.5 3 3.5 4 

mass (A/©) 


Figure 5. Calculated ratios for the 21 giants after 

FDUP (solid l ine) and SDUP (dashed line) versus stellar mass . 
Cal culation by l B oothroyd fc S ackmannl lll999l) . [Abia et ^ ll2012l) 
and iKarakas &; Lattanziol l|20lJ) are also shown for comparison. 
Observations are indicated as data points. The errors on the 
masses are as listed in Table [2] 


extra mixing or overshooting at the bottom of the convective 
envelope. In order to constrain this overshooting, a better 
understanding is required for the role of first and second 
dredge up for the whole range of masses under consideration. 

It is known that after the star leaves the main sequence, 
FDUP alters significantly the surface abundances as it en¬ 
riches the envelope with “^He, and and reduces 

its and abundance. Every star that evolves 

to the AGB experiences the FDUP episode, and starts its 
early AGB phase with a sharp composition discontinuity at 
the point of maximum penetration of the FDUP. The H- 
burning shell in low-mass stars represents an entropy bar- 
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Table 3. (in solar units) and the surface number abundances after FDUP and SDUP (if any, listed in parenthesis): a: surface 

abundances with standard mixing (present work), 6 : TsujiOS, 3 : Harris and Lambert (1984b), h: H85: Harris et al. (1985), v. Smith and 
Lambert (1990), j'. Decin et al. (2003), fc: Maillard (1974), Z: Hinkle et al. (1976). 


Object 

Mhgs 


16o/17o<> 

160/170 

120/130^ 

12c/13c6 

12c/13c 

14n/15N“ 

5 Vir 

1.19 

2370 

> 2500 


32 

12.3 ± 1.2 

16 ±4* 

551 

a. Tau 

1.2 

2370 

> 1000 

^'-'^+130 

30 

10.6 ± 1.0 

10 ± 2hJ 

551 





K9i:;~125g 

^^^+250 



9± 19 


RRUMi 

1.15 

2006 

> 2000 


30 

10.0 ±0.8 


598 

RZ Ari 

1.14 

2006 

607 ± 48 


30 

7.9 ± 0.8 


598 

5 Oph 

1.39 

1500 

387 ± 68 


29 

11.1 ±0.9 


690 

V Vir 

1.39 

1500 

> 2000 


29 

8.7 ± 1.3 

12 ± 2* 

690 

Eri 

1.73 

461 

687 ± 14 


24 

12.4 ±0.3 


963 

10 Dra 

1.74 

461 

151 ± 11 


24 

14.8 ± 1.6 

12 ± 3“ 

963 

/3 Peg 

1.7 

462 

> 2500 

1U0U_|_5QQ 

24 

7.7 ±0.5 

8 ±2* 

963 





3s 1009 



5 ±39 


3 O 3 Her 

1.65 

462 

211 ±42 


24 

12.5 ± 1.1 

10±2 * 

963 

a Lib 

1.9 

301 

> 1500 


22 

7.5 ± 0.3 


1126 

R Lyr 

1.8 

301 

368 ± 44 


22 

6.4 ±0.3 


1126 

fi Gem 

1.9 

301 

798 ± 73 

09c —75 g 

'J^'J+150 

22 

10.5 ± 1.2 

13 ± 2* 

1126 

3: 1009 

OP Her 

1.9 

246 

329 ± 31 


22 

11.3 ± 1.2 


1103 

p Per 

2.4 

234 

> 1000 


22 

9.7 ± 1.0 

15 ± 2' 

1307 

a. Get 

2.96 

318(317) 

586 ± 47 


21 (21) 

11.1 ±0.8 

10 ± 29 

1460(1501) 

A Aqr 

2.96 

318(317) 

> 1000 


21 (21) 

7.9 ± 1.4 


1460(1501) 

XY Lyr 

3.0 

424(395) 

223 ± 16 


21 (20) 

15 ±0.4 


1528(1610) 

52 Lyr 

4.5 

400(378) 

465 ± 41 


21 (20) 

16.2 ± 1.5 


1503(1629) 

ol Her 

4.0 

424(402) 

102 ±8 

i8o;?°« 

21 (20) 

11.1 ±0.7 

17 ±494 

1538/1642 

200 +^®'' 

Rri 4509 

4e;n~50fc 

40U^50 

BS 6861 

5.35 

424(402) 

> 1000 


21 (20) 

48.5 ±2.9 


1538(1642) 

Initial 


2620 



90 



270 


rier that prevents any deeper penetration of the convective 
envelope, and thus no further change in the surface abun¬ 
dances takes place. However, for solar metallicity stars of 
masses above 3-4M0, the gravitational energy release due 
to the contracting core and the increased energy flux from 
the burning He-shell lead to an expansion so that the H-shell 
is pushed outwards in mass to low temperatures causing a 
temporary extinction in the H-burning shell. This situation, 
coupled to the drop in the temperatu re of the expanding lay - 
ers and the increase in the opacity (llben fc R.enzinilll98,3l l. 
causes convection to deepen in a second dredge up (SDUP) 


event and mixes out the composition discontinuity left over 
by the FDUP. Therefore, SDUP introduces further changes 
in the surface abundances of stars of masses > 4Mq. 

Fig. [6] shows the maximum penetration of the convec¬ 
tive envelope at FDUP and SDUP as a function of the initial 
stellar mass. It is clear from the figure that for M< (4-5)Mq, 
SDUP doesn’t penetrate deeper than the FDUP, and thus, it 
does not introduce signihcant c hange to the surface compos i- 
tion. This is in agreement with iKarakas fc Lattanzi^ ll2014h . 
for their solar metallicity case. Their Fig. 7 exhibits similar 
general features and also indicates a deepest penetration of 
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M(Mq) 

Figure 6. Fractional mass reached by the convective envelope 
at FDUP (dashed line) and SDUP (solid line). SDUP is deeper 
that FDUP only in models of stellar masses SMq. The mass 
on the y-axis is given as the mass reached by the dredge up over 
the total initial stellar mass, in solar units. 



mass (A/q) 

Figure 7. The ratios obtained in the standard case 

and in the case of envelope overshooting with /=0.125. Several 
observations are shown for comparison. 


FDUP in « 2 . 5 M 0 stars. For stars of M> 6 M 0 , the expan¬ 
sion is strong and convection extends deeper inwards. Since 
low mass stars do not experience SDUP, this implies that 
their observed surface abundances on the AGB phase are 
actually “inherited” from the RGB phase. Therefore, one 
way to account for the discrepancy between calculated and 
observationally inferred values of in these stars is 

to invoke extra mixing below the convective envelope on the 
RGB phase. 


3.2.2 Predictions of surface abundances with extra mixing 


Overshooting is applied as outlined in Section 2.1 and we 
Hnd that / = 0.125 provides the best estimation for its ef- 
hciency in low mass stars. Fig. [7]shows the ratios 

of our star sample in the standard case and with overshoot¬ 
ing. The low observational in low mass stars can 

now be fitted quite well within the error bars on the stellar 
masses. It is noted how overshooting has a minor effect in 
stars of M 3 M 0 . This is expected due to the shape of the 
abundance profile in these stars as discussed earlier in 
connection with the shape of the prohle in these stars. 
The calculated and observationally inferred values for these 
stars are generally in a good agreement. 

Fig. [S] shows the predicted carbon and nitrogen abun¬ 
dances in the standard case and with overshooting. Observa¬ 
tions are also shown for comparison. A well-known problem 
arises in explaining the surface carbon abundances in stars of 
2 M 0 , where our extra mixing treatment cannot explain 
the low carbon observ ed in these giants. Calculation s us- 
ing standard mixing bviBoothrovd fc SackmanrJ lll999l) and 
iKarakas fc Lattanzid (12014 ) show a similar problem. A dif¬ 


ferent non-standard mixing process seems to be required in 
low-mass stars (< 2 M 0 ) which can reduce the abundance 
of *^^0 by mixing it to the hotter regions of the H-burning 
shell, allowing for some nuclear processing. Such mixing pro¬ 


cess is linked to the long evolutionary time of the low-mass 
st ars while ascending th e RGB . As discussed in some details 
in IKarakas fc Lattanzid (l2014f) . the so-called “thermohaline 
mixing” or “salt-hnger instability” seems to be one possible 
description of this extra mixing in low-mass red giant stars. 
This is still an unsettled problem in stellar modeling. Nitro¬ 
gen abundances, on the other hand, show a scatter in the 
observed data. This will be discussed in connection with the 
i^N/isn ratios in Section 3.3.2. In general, Fig. [SJa shows 
that the observations pose no serious contradiction with the 
values predicted by the stellar models which, on average, 
reasonably reproduce the observationally inferred data. 


3.3 Effect of modified nuclear reaction rates 

3.3.1 The OO ratio 


In order to investigate the effect of the "O production 
and destruction reaction rates on the ratios, 

four different evaluations of the proton-capture reactions 
i® 0 (p, 7 )^'^F, i^O(p, 7 )^®F and ^^0 (p, q)^'^N are u s ed. In 

(|2014|) 

^ 20101 ) 


particul ar, we use the comp ilations by lSergi_et_^ 


(SE14), ISallaska et al.l (l2013f) (SA13). Illiadis et al 


(ILIO) and iGhafa etaP (l2007l) (CH07). The ra¬ 


tios obtained are shown in Fig. [§] for both cases: standard 
mixing (Fig. |9^) and envelope overshooting with / = 0.125 
(Fig. [^). The four sets of rates give very similar 
values in the considered mass range. None provides a reason¬ 
able agreement between model predictions and observations 
unless overshooting is included. This consistency shows that 
the existing discrepancy cannot be removed without invok¬ 
ing deeper mixing. Fig.l^j shows a better fit of the observed 
data within the error bars in the low mass stars. 

The effect of the reaction rates uncertainties is also 
worth exploring. Since the SA13 compilation is based on 
a Monte Carlo simulation and the rates have s t atistically 
well-defined uncertainties (jLongland et al.l l201fll : iLongland 
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Figure 8. The predicted surface abundances of (a) carbon and (b) nitrogen with observations as shown. The initial solar values of carbon 
and nitrogen are 8.5601 and 8.0499, respectively and are indicated by a star symbol. 




Figure 9. The ratios calculated for the 21 giants, with four different sets of compilations of relevant reaction rates (see text 

for details). Left panel shows the case with standard mixing and the right panel is that with overshooting (/ = 0.125). Observational 
data are also included. 


|2012| : liTiadis et al.ll2014h . the ratios are calculated 

using the recommended, high and low rates, where the rate 
boundaries correspond to a 95% coverage probability. The 
rates uncertainty has a very minor effect on the tracks during 
H-shell burning but none along the RGB, and hence it does 
not affect our mass determination. However, the 
ratios show a larger sensitivity to these uncertainties, as 
shown in Fig. HOI It is found that the observational data can 
be better explained when the whole range of uncertainty on 
the involved reaction rates is considered since the discrep¬ 
ancy between predictions and observations becomes less pro¬ 
nounced. On the other hand, the ratios are found 

to be immune to the uncertainties on the considered rates. 
The difference in obtained with the high rates and 

low rates does not exceed 3%. 


3.3.2 The ratio 

The surface ratio is worth considering in connec- 

tion with the late st evaluation of the ^'’N(p, 7)^®0 rate 
llMarta et aLllioiIh . This ratio is calculated for the mass 
range under consideration after FDUP and SDUP. It is 
fo und tha t our v alues are higher by ~ 20% relative to those 
bv lEl EidI (|l994h in stars of masses SM©. T his is expected 
and due to the NACRE ^'*N(p, 7)^®0 rate jAngulo et al.l 
Il999l) used in that work, which is almost a factor of 2 
higher than the revised rate at ste llar temperatures (see 
Halabi. El Eid fc Champagnel ll2012h for the expression of 
this rate and explicit discussion). However, our is 

about 12% lower in the low-mass stars. 

It is quite unfeasible to verify our findings since 
ratios are difficult to measure in RGB stars from CN lines 
because they are too weak even in very high resolution spec¬ 
tra. In fact, the ratio is difficult to measure di- 
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Figure 10. The SA13 set of rates in the standard case. Shown 
are the results with the recommended set of rates (solid line), the 
high rates (dotted) and low rates (dashed). Observational data 
are the same as in Fig. |9] 


rectly for molecules containing carbon llRoueff et al.l[2015l) . 
On the other hand, stars in the early AGB phase are O- 
rich and then their spectra would be dominated by oxide 
molecules and not by C-bearing molecules. Oxide molecules 
like NO cannot be observed in the visual or near-IR region 
in these stars (Carlos Abia, private communication). Ob¬ 
servations in radio-wavelengths can in principle detect some 
N-bearing molecules from which may be derived, 

but in this case one would probably be looking at the cir- 
cumstellar envelope of the star not the photosphere. Cir- 
cumstellar N ratios might be affected by the incoming UV 
radiation from the ISM triggering non-kinetic equilibrium 
chemistry and thus, might not rep resent the stellar photo- 
spheric ratios (iHedrosa et al.ll2012l) . Theoretical predictions 
of the ratios, particularly the changes induced by 

the revision of the ^^N(p, 7)^®0 reaction rate may benefit 
from future observations or one may resort to the chemical 
analysis of pre-solar grains originating from the envelopes of 
late AGB stars. 


4 CONCLUSIONS 

A sample of observed RGB and early AGB stars was con¬ 
sidered and their masses were obtained using extended evo¬ 
lutionary tracks. This was done without the need of extrap¬ 
olating the evolutionary tracks at lower effective tempera¬ 
tures as it has been done in the work by TsujiOS. The present 
evolutionary tracks include the effect of mass-loss, which be¬ 
comes important during the red giant phases, especially for 
the higher stellar masses under consideration. This investi¬ 
gation includes an analysis of the physical conditions in the 
overshooting region and the effect of this overshooting on 
the abundance yields. We find that overshooting is needed 
to reconcile observational oxygen abundances with model 
predictions particularly in low mass red giants. Although 


TsujiOS has interpreted the discrepancies between the pred¬ 
icated and observational abundances as an effect of extra 
mixing, he did not provide models including this extra mix¬ 
ing to see how the discrepancies may be understood. We 
showed that extra mixing based on the mixing of chemical 
elements by diffusion explains reasonably the observation- 
ally inferred oxygen isotopic ratios. It is important to realize 
the challenges facing such observations and the uncertain¬ 
ties involved in the available data. In this regard, the spread 
in the observational data in the low and high mass ranges 
was discussed in connection with the inherent difficulties in 
analyzing the spectra of these relatively cool stars and the 
uncertainties involved in measuring faint lines like 
Our overshooting treatment cannot, however, explain the 
low surface carbon abundances in low mass stars. Another 
mixing mechanism seems to be required during the long evo¬ 
lutionary time needed for the low-mass stars to ascend the 
red giant branch. 

Furthermore, the present calculations were carried out 
using recent determinations of proton-capture rates which 
have reliable statistical error bars. This allows us to draw 
conclusions on the uncertainties involving GNO surface 
abundances. In particular, the effect of recent evaluations 
of the reaction rates on the production and destruction of 

was explored. The experimentally suggested uncertainty 
of these rates provides a better fit of the observed 

in low mass stars yet does not exclude the need to invoke 
overshooting. Additional mixing beyond the convective en¬ 
velope as determined by the Schwarzschild criterion is found 
to be necessary to better explain the observational 
surface abundances, especially in low mass stars. Moreover, 
the effect of ^^N(p, 7)^®0 rate on the ratios was 

studied in the considered mass range. 

As a final remark, our approach in the present study 
was to consider a sample of red giant stars to see how far 
standard calculation predictions agree with observations. A 
comprehensive comparison between stellar models and ob¬ 
servations based on the analysis of the effect of extra mixing 
was presented and linked to the uncertainties in key nu¬ 
clear reaction rates affecting the GNO abundances in red 
giants. Comparing theoretical predictions of stellar models 
to available observations is required in order to constrain 
parametrized approaches in determining the efficiency and 
extension of mixing at convective boundaries. Future multi¬ 
dimensional simulations of convection may introduce an im¬ 
proved local description of this mixing and provide insight 
towards a better understanding of the physical processes in¬ 
volved. 
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the manuscript. The authors also acknowledge the anony¬ 
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